arXiv:1507.01594vl [hep-ph] 6Jul2015 


MCTP-15-14 


Status and prospects of precision analyses with e^e — W~^W 

James D. Wells and Zhengkang Zhang 
Michigan Center for Theoretieal Physies, Department of Physics, 

University of Michigan, Ann Arbor, MI 4^109, USA 

(Dated: July 8, 2015) 

e^e“ —7- W^W~ is an integral part of the global precision analysis program which 
is becoming more relevant after the discovery of the Higgs boson. We analyze the 
current situation of precision calculations of inclusive e'*“e“ —)■ W^W~ observables, 
and study the prospects of incorporating them into the framework of global pre¬ 
cision electroweak analyses in light of per-mil-level cross section measurements at 
proposed future facilities. We present expansion formulas for the observables, mak¬ 
ing the dependence on the inputs clear. Also, the calculation of new physics effects is 
demonstrated in the effective held theory framework for universal theories. We go be¬ 
yond the triple-gauge-couplings parametrization, and illustrate the complementarity 
of e+e” —)■ W^W~ and other precision data at the observables level. 
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I. INTRODUCTION 

The recent discovery of the Higgs boson BEI has initiated intense interest in precision 
stndies of its properties. There have been various attempts to connect Higgs phenomenology 
to previous precision electroweak measurements, in order to see what windows for new 
physics are still open and can be probed by precision Higgs measurements |3HI3]- There are 
two facts behind this trend. First, deviations from the Standard Model (SM) in the Higgs 
and electroweak sectors are correlated if the recently discovered Higgs boson is part of an 
SU{2)l doublet. To fully exploit the power of precision analyses we should take advantage 
of all the strong constraints on new physics from precision electroweak studies. Second, 
it has been realized that in the effective field theory (EFT) framework, which is the only 
consistent model-independent approach to new physics with minimal assumptions [HHTB], 
it is important to use a complete operator basis [ini [la IIM9]. Therefore, previous EFT 
calculations of precision electroweak observables should be recast into a uniform framework 
for a global analysis. In this context, a crucial role is played by e+e“ —)■ W^W~. Despite 
the lower experimental precision compared with the most precisely measured observables, 
—)■ W^W~ offers a unique window to test the SM and probe new physics effects due 
to its sensitivity to triple-gauge couplings (TGCs) which are difficult to access otherwise. 

On the other hand, precision studies of H^-pair production are important in their own 
right. Such studies have been carried out for collisions at LEP2 pOl - EH] . for pp collisions 
at the Tevatron [23], and more recently for pp collisions at the LHC [221 EH]- The clean 
experimental environment at LEP2 allowed e’''e“ —)■ W^W~ cross sections to be measured 
at the ~ 1% level up to ^/s = 207 GeV, and agreements were found with SM predictions 
with a similar precision [271138] . Historically, analyses have been done in the language of 
TGC parameters [39]. In this context, LEP2 data still provide the most stringent constraints 
on anomalous TGCs, with LHC measurements just starting to become competitiv^ Future 
high-luminosity colliders, operating at center-of-mass energies from W^W~ threshold 
up to 1 TeV or even beyond, will enable per-mil-level cross section measurements in a wide 
range of y/s and push the precision frontier on —)■ W^W~ studies much further [35L 

139] . Accordingly, progress has been made after LEP2 on precision calculations toward the 
per-mil precision goal [501 El]. The promising experimental progress calls for reassessment 


^ The tPy and WZ channels are also studied at the LHC to put limits on anomalous TGCs [3(11173] . 
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of the role of e’''e —)■ W~^W in the precision program, both as a consistency test of the 

SM and an indirect probe of new physics. 


In this paper, we revisit the calculation of several —)■ W~^W~ observables, both in 

the SM and in the presence of new physics, in the expansion framework of Ref. [52] . From 
the point of view of testing the SM, at present our best knowledge of the compatibility of 
the electroweak SM with data comes from global analyses of Z-pole observables and miy, 
all of which have been very well-measured [53] [M] and precisely-calculated (55] |56]; see e.g. 
Refs. [HI 1561160] for recent global hts. It is interesting to ask whether improved measurements 
of —)■ W~^W~ will play a complementary role in such global analyses (aside from a 

better determination of my/). To answer this question we need to gather information on 
how the —)■ W^W~ observables depend on the inputs of the calculations, which has 

not received much attention in the past. The expansion formulas we present will make the 
answer transparent. 


To demonstrate the role of e’''e“ —)■ W^W~ in probing new physics, we adopt the EFT 
approach. We will focus on an important class of new physics scenarios, the so-called “uni¬ 
versal theories” [6l], for illustration. The results will be presented in a way that allows other 
precision constraints to be easily incorporated. Many previous studies used the reported ex¬ 
perimental values for the TGCs to constrain the EFT parameter space, but there is a caveat 
related to the assumptions made when extracting the TGCs [ 13 ]. In this regard. Ref. [ 13 ] 
refers to the TGCs as “constructed observables”, and points out that extreme care must be 
taken when relating constructed observables to EFT parameters. Here we take a different 
approach by working with well-dehned physical observables which are free from such sub¬ 
tleties. A similar analysis has been done in Ref. [62] without discussing the interplay with 
other precision measurements. 


We will begin in Sec. [^by dehning several e+e —)■ W~^W observables and reviewing 


their calculations in the SM. Sections III and IV are devoted to the SM and new physics 


aspects of precision analyses mentioned above, respectively. In Sec. |V|we conclude. 
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II. OBSERVABLES AND STANDARD MODEL CALCULATIONS 


At leading order (LO), the process —>■ W~^W~ is calculated in the SM from the 

,s-cha nn el Z /7 exchange and t-channel neutrino exchange diagrams, known as the CC03 
diagrams | 2 ll | 22 ]. Some LO results that will be used later are collected in Appendix [A} 
However, the notion of —?• W~^W~ observables” is not well-dehned once we go beyond 

LO to include hnite lU^ width effects and radiative corrections. Strictly speaking, what 
we refer to as “e+e“ —>■ W~^W~ observables” should be understood as a shorthand for 
“e+e“ —>■ 4/ observables with 4/ compatible with intermediate W~^W~”. Even at tree 
level, there are diagrams involving only single intermediate W~^ or W~ that contribute to 
e’''e“ —)■ 4/. They are known as CCll diagrams that are not in the CC03 class [211 [ 22 ], ^-nd 
interfere with the CC03 diagrams with intermediate W~^W~. We will include all the CCll 


diagrams in the calculation of SM predictions in Sec. |III[ as opposed to the calculations 
adopted by LEP2 analyses where only CC03 diagrams are included [23]. One should keep 
in mind that with hnite width, only the sum of the complete set of CCll diagrams is 
gauge invariant. 

The —)■ W^W~ observables we consider are the polarized total cross sections ctl, aji, 

dehned as the cross sections with left- or right-handed incoming electron and unpolarized 
incoming positron. From these, the unpolarized total cross section 


(T = -{(Tl + ctr), 


and the left-right asymmetries 


_ ai - aR _ 
^LR = - — t 


-, ^LR = t 

a 


4 - 

4iLR — - 


( 1 ) 


( 2 ) 


o-L + o-R 0- a 

can be derived. Aur may be more convenient than Alr because the latter is very close to 
1. We will focus on these inclusive observables without applying kinematic cut^ Other 
observables can be extracted from differential distributions, which can presumably be dis¬ 
torted by new physics effects. For example, the forward-backward asymmetry dehned for 
the outgoing W~ is often considered; see e.g. [62] for an EFT study of diherential cross 
sections at LO. However, once nonfactorizable radiative corrections and oh-shell ehects are 


^ The numerical difference between our reference values for cr in Sec. 
mostly due to the separation cuts applied in the latter. 


Ill 


and the results in Refs. 


IS 
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taken into account, experimental subtleties arise related to the kinematic reconstruction of 
from the 4/ hnal states, which should be carefully studied and is beyond the scope of 
the present paper. 

The state-of-the-art calculations of e’''e“ —?• W~^W~ cross sections in the LEP2 era in¬ 
corporated O {a) radiative corrections in the double-pole approximation (DPA), and were 
implemented in dedicated programs RacoonWW [27H33] and YFSWW3 [541 - I5B] . Later 
the complete O {a) radiative corrections were calculated for the four-fermion hnal states 
i/rT^udfi~Ufj, and udsc isni EH- However, the latter calculation is not yet available 
as public codes. So for the purpose of illustration we will present the results as calculated 
within the DPA implemented in the program RacoonWW [33]. The unpolarized total 
cross section in the DPA agrees with the complete O (a) result within 0.3% for ^/s from 
200 GeV to 500 GeV [SU]. To achieve better precision suitable for studies at future colliders, 
the results presented here are expected to be updated once more up-to-date codes become 
available. We also remark that to achieve better-than-per-mil accuracy, even the complete 
O (a) calculation needs to be supplemented by higher-order Goulomb corrections near the 
threshold, higher-order Sudakov logarithms for a/s > 500 GeV, and QGD effects [4611501163] . 

We will consider two benchmark center-of-mass energies ^/s = 200 GeV and 500 GeV, 
where the DPA works reasonably well. Among all the four-fermion hnal states compati¬ 
ble with intermediate W~^W~, we will focus on udfi~9^ for illustratio]:]|^ This channel is 
representative of the mixed leptonic and hadronic decay channels from W^W~, which are 
expected to have high selection efficiency and a low background. Separate calculations may 
be needed if one is interested in other channels due to additional diagrammatic contributions 
and the necessary inclusion of hnite-electron-mass effects in the case of hnal-state e^. 


III. PARAMETRIC DEPENDENCE AND UNCERTAINTIES 

We will present the results of the SM calculations in the form of expansion formulas 
as in Ref. [52], which make clear the parametric dependence and uncertainties. To briehy 
review the formalism, we note that the SM predicts each observable as a function of 
^ Theory calculations are conventionally formulated in terms of partonic final states, though experimentally 
jets instead of quarks are observed. In this sense what we call “observables” are not yet experimentally 
observed quantities, but are directly related to the latter when we sum inclusive quark contributions that 
form jets. 
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the Lagrangian parameters. A more convenient way to arrange the calculation, which is 
commonly adopted in precision electroweak analyses, is to eliminate the input Lagrangian 
parameters in favor of the same number of very well-measured observables {O-,} (as in 
Ref. [52], we use primed indices for input observables and unprimed indices for output 
observables). Then the analysis involves only physical observables, so that all results have 
unambiguous meanings, and comparison with experiment is straightforward. The calculation 
is further simplihed by expanding the theory (SM) predictions about some reference point 
and keeping terms up to linear order: 


= Of^ = Of 


E 


ddfM 


f do. 


where 


8.,=8’^f 


Of 


d, - df) = df (i + Y, ), (3) 


( 4 ) 


is the fractional shift of the input observable O, with respect to its reference value, and 


Oa' — 


df ddf^ 


of do. 


( 5 ) 


d/=oT 


represents the resulting fractional shift of the output observable O^ calculated in the SM. 
The superscript “SM” in Eq. (Q indicates that the shift can be associated with adjusting 
the SM Lagrangian parameters, which should be distinguished from corrections due to new 


physics; cf. Eq. (22). These expansion coefficients characterize the parametric dependence 
of the calculated observables on the input observables, as long as the expansion up to hrst 
order is adequate. This is the case for most practical purposes now that the mass of the 
Higgs boson is known to subpercent level, and so all have to be much smaller than 

unity. 


The standard set of input observables commonly adopted in precision electroweak analy¬ 
ses consists of the masses of the Z boson, the top quark, and the Higgs boson 
the Fermi constant Gp, and the couplings a{mz),(ys{mz). In the case of e+e“ —)• W~^W~, 
however, it is more convenient to use the W boson mass mw in place of a{mz) as a cal- 
culational input. Thus, we hrst extract the expansion coefficients with respect to the input 
observables 


{Of} = 


(6) 









and then transform the resnlts into the standard basis 


= {'mz,GF,a{mz),'mt,as{mz),mH}- 


( 7 ) 


The SM predictions for the observables take the following forms in the two basis: 



or 

or 


+ Ci^asS^^Ctsimz) + 

1 + + di^CpS^^Gp + diA^^cd'rnz) 

+di,mj^^mt + di^aj^^asimz) + 


The transformation from the Cj_j/ coefficients to the diy coefficients 


di^a ^i,m\^rdm\iY ,Ofi 

di^i' Cj,i' T dn^^y ji' for Oj/ G \tTI z •! G p (y,g{TTl z) •!'Tflpi^ 

can be derived similarly as in Section 3.4 of Ref. |52j . 

We adopt the following reference valnes for the inpnt observables: 

mf = 91.1876 GeV, 

Gf = 1.1663787 x 10"^ GeV^ 
mr = 80.3614 GeV, 
mr = 174.17 GeV, 

= 0.1185, 
mr = 125.9 GeV. 


These, according to the formulas presented in |52], correspond to 

a{mzA^ = 7.81861 x 10“^ = 1/127.900. 


( 8 ) 

(9) 


( 10 ) 

( 11 ) 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 


(18) 


The hnal results of the SM predictions for the hve e’''e“ —?• W^W~ observables are the 
expansion formulas Eqs. ([^ and ([^, with the reference values and expansion coefficients 
listed in Tables and respectively, for the benchmark center-of-mass energies a/s = 
200 GeV, 500 GeV. The error bars quoted contain (possibly overestimated) Monte Garlo 
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vG/GeV 

o. 



G,Gf 

^i,mw 


^i,mH 


cTL/fb 

1229.8(5) 

0.123(15) 

1.957(15) 

1.429(15) 

0.0038(29) -0.0141(29) 

0.0006(29) 


CTi^/fb 

13.593(5) 

17.48(8) 

1.981(13) 

-17.03(8) 

0.1364(27) -0.0138(27) 

-0.041(5) 

200 

tr/fb 

621.70(25) 

0.312(15) 

1.958(15) 

1.227(15) 

0.0052(29) -0.0141(29) 

0.0002(29) 


Alr 

0.978136(12) 

-0.3835(18) 

-0.0005(4) 

0.4078(18) 

-0.00293(9) -0.00001(9) 0.00092(14) 


Alr 

0.021864(12) 

17.17(8) 

0.024(20) 

-18.26(8) 

0.131(4) 0.000(4) 

-0.041(6) 


CTL/fb 

553.48(22) 

0.341(18) 

2.022(18) 

2.936(19) 

0.0023(29) -0.0120(29) 

0.0005(29) 


CTfi/fb 

3.4687(13) 

14.93(7) 

2.098(12) 

-11.04(7) 

0.1710(25) -0.0104(25) 

0.0042(25) 

500 

cr/fb 

278.48(11) 

0.432(18) 

2.022(18) 

2.849(19) 

0.0034(28) -0.0120(28) 

0.0005(28) 


Alr 

0.987544(7) 

-0.1826(10) 

-0.00096(28) 0.1751(10) 

-0.00211(5) -0.00002(5) 

-0.00005(5) 


Alr 

0.012456(7) 

14.49(8) 

0.076(22) 

-13.89(8) 

0.168(4) 0.002(4) 

0.004(4) 


TABLE I. Reference values and expansion coefficients for the e+e“ —)• W~^W~ observables (for 
final state) with respect to the input observables {mz,GF,rnw,'nT't:Cts{mz),mH}, to be 
used in Eq. Q. 


vG/GeV 

o. 

Off 

di^raz di^Gp 

^Z,Q( 

^i,mt ^i^OLs di,mH 


CTL/fb 

1229.8(5) 

2.163(26) 2.272(16) 

-0.3077(33) 

0.0227(29) -0.0154(29) -0.0005(29) 


CTfi/fb 

13.593(5) 

-6.84(14) -1.767(22) 

3.668(17) 

-0.0893(29) 0.0026(27) -0.028(5) 

200 

cr/fb 

621.70(25) 

2.064(26) 2.228(15) 

-0.2643(32) 

0.0215(29) -0.0152(29) -0.0008(29) 


Alr 

0.978136(12) 

0.1988(32) 0.0892(6) 

-0.0878(4) 

0.00247(9) -0.00040(9) 0.00061(14) 


Alr 

0.021864(12) 

-8.90(14) -3.995(27) 

3.933(17) 

-0.111(4) 0.018(4) -0.027(6) 


CTL/fb 

553.48(22) 

4.534(32) 2.668(19) 

-0.633(4) 

0.0412(29) -0.0148(29) -0.0017(29) 


CTfi/fb 

3.4687(13) 

-0.85(13) -0.333(20) 

2.379(16) 

0.0246(27) 0.0002(25) 0.0127(25) 

500 

cr/fb 

278.48(11) 

4.501(32) 2.649(19) 

-0.614(4) 

0.0411(28) -0.0148(28) -0.0017(28) 


Alr 

0.987544(7) 

0.0674(17) 0.03757(35) 

-0.03771(21) 

0.00021(5) -0.00019(5) -0.00018(5) 


Alr 

0.012456(7) 

-5.35(13) -2.982(28) 

2.993(16) 

-0.016(4) 0.015(4) 0.014(4) 


TABLE 11. Reference values and expansion coefficients for the e'’'e“ —)• W^W~ observables (for 
final state) with respect to the input observables {m^, Gi?, a(m 2 ), as(m 2 ), m/f}, to 

be used in Eq. 
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errors only, while truncation errors from numerical differentiation are expected to be smaller. 
Further technical details of the calculation can be found in Appendix 

To get an idea of the size of O (a) radiative corrections, as well as the dependence on 
yjs beyond the two benchmark choices, we show in Fig. [^the comparison of the numerical 
results in Table with LO results as functions of y/s for the two observables and a^. 
The latter can be easily calculated analytically; see Appendix 0 Only di^rnz, di^Gp, 
di^a are presented, since the other three expansion coefficients vanish at LO. Fig. shows 
that in most cases, the O (a) corrections are at or below O (10%) level. We note that while 
detailed discussions on the size of O (a) corrections for can be found in the papers 
where these corrections are calculated, the plots for the expansion coefficients di^i' are new. 
The latter provide complementary information because both and di^ii are needed in 
precision analyses if one is not restricted to fixed values of the input observables. Also, di^i' 
are essential in a consistent calculation of new physics effects, as we will see in the next 
section. 

The results in Tables |I| and |IT| reflect the accuracy implemented in Racoon WW, namely 
up to O (a) (with respect to LO) and within the DPA. They are expected to be updated 
in the future. However, even at present, these finite-accuracy expansion formulas are useful 
for the purpose of having a picture of parametric dependence and an estimate of para¬ 
metric uncertainties, i.e. uncertainties from the input parameters (observables). With the 
experimental uncertainties of the input observables taken from Ref. [S2], 


Airiz = 2.1 MeV, AGp = 6 x 10 GeV Aa{mz) = 8.6 x 10 

Amt = 0.87 GeV, Aas{mz) = 0.0006, Amn = 0.4 GeV, (19) 


the fractional parametric uncertainty in from input observable is easily obtained by 



I di^i 


AO, 

df 


= At,' ■ 10 


-4 


( 20 ) 


The results for are listed in Table III for two representative observables aR and a. 


The parametric dependence and uncertainties for —)■ W~^W~ observables are usually 

not discussed in the literature, but they provide important information if we put these 
observables into the broader context of precision electroweak analyses. In particular, a global 
fit is dominated by observables for which experimental and theoretical uncertainties 
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FIG. 1. Comparison of the RacoonWW calculations (red data points) of and expansion 
coefficients with LO results (continuous curves) as functions of ^/s. 
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V^/GeV 

o. 

^i,mz 

^i,GF 


^i,mt 

^i,as 



CTR 

1.57 

0.01 

4.07 

4.46(15) 

0.13(14) 

0.89(17) 

200 

a 

0.48 

0.01 

0.29 

1.07(14) 

0.77(15) 

0.02(9) 


ctr 

0.19 

0.00 

2.64 

1.23(14) 

0.01(13) 

0.40(8) 

500 

a 

1.04 

0.01 

0.68 

2.05(14) 

0.75(14) 

0.05(9) 


TABLE III. Fractional parametric uncertainties from each input observable, in units of 10“'^. For 
example, at ^/s = 200 GeV, ^ajo from raz is 0.48 x lO”'^. These parametric uncertainties are 
negligible compared with experimental and theoretical uncertainties. 

are not much larger than parametric uncertainties, e.g. the effective electroweak mixing 
angle sin^ In this regard, we hnd that the projected per-mil-level measurements and 
calculations of the e+e“ —)■ W^W~ total cross section a are still not good enough compared 
with the very small parametric uncertainties. The latter are seen from Table Ell to be 
at the 10“^ level at present, and will be further reduced in the future with more precise 
measurements of the input observables. The parametric uncertainties are larger for cx/j, 
but are still expected to be much smaller than the experimental errors, given the limited 
statistics for an since aji a. We thus conclude that the —>■ W^W~ observables 

are likely to remain of peripheral importance in a global precision electroweak analysis as 
a compatibility test of the SM. Nevertheless, —)■ W~^W~ can be a uniquely powerful 

probe of some new physics scenarios, as we will describe below. For this purpose, the results 
in this section indicate that new physics effects above the O (10“^) level can be studied 
without letting the SM input observables float even in a global analysis. 


IV. PROBING UNIVERSAL THEORIES WITH e+e" ^ W+W- 

A. Universal theories and operator bases 

In the presence of new physics, an additional term is added to the right-hand-side 
of Eq. (|^, which is dehned as the fractional shift in Of" due to new physics calculated in 
terms of the Lagrangian parameters. It is assumed that |.^i| -Cl so that an expansion up to 
linear order is still adequate. To conform with the very precise measurements of the input 
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observables, we adjust the SM Lagrangian parameters to keep 

- df 


th _ 


f = 






( 21 ) 


small. Then the output observables are calculated as follows: 


0 “ = of b + ^ + 5.) = of M + ^ d„,sdf + s"'’ob, (22) 


where 


^ ^ ^ di^rnz^mz difi^^Qp di^a^o, 


(23) 


In Eq. (22), the shift in with respect to is consistently decomposed into the shift 


due to tuning the values of the input observables and the shift due to new physics. The 
latter is in turn decomposed into a direct contribution and an indirect contribution from 
undoing the shifts in the input observables. 


We will apply this formalism to a popular class of new physics scenarios, the “universal 
theories”, and investigate their effects on e^e“ —)■ W^W~. Universal theories are dehned by 
the assumption that new vector states, if there are any, couple to SM fermions only via the 
SU (2) L X U{1)y currents [HI]. In other words, it is assumed that the only gauge interactions 
of the SM fermions, apart from QCD, have the form 

A£ = 9ll'^“((ryi + g7'‘y,)+g'S„ F,/V/. (24) 

f&{l,e,q,u4} 

where IU“, need not coincide with the SM gauge helds. Simple examples that qualify 
as universal theories include W, Z' models, where IU“, are mixtures of SM and new 
vector bosons. More complicated new physics models that are well-motivated (e.g. little 
Higgs models and some extra-dimension models) can also be cast into this form |61j . 

If we further assume that the scale of new physics A is somewhat higher than ^/s (which 
is well-motivated for y/s ^ 500 GeV given the non-observation of particles beyond the 
SM so far, though the situation might change), their effects can be parametrized model- 
independently by an EFT, 

C = Csm + Y^ ^Oi + O where ~ O . (25) 

By organizing the EFT as an expansion in n^/A^, we have assumed the recently-discovered 
Higgs boson is part of an SU{2)l doublet. Relaxing this assumption leads to the nonlinear 
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version of the EFT, which is more appropriately organized as a loop expansion iumdsi 
EHES]- We will calculate at LO, namely O terms arising from the interference 

between tree-level diagrams with one dimension-six operator insertion and LO diagrams in 
the SM. In particular, this means CP-odd operators do not enter. Also, in the limit —?• 0, 
universal theories can modify W~^W~ production but not W~^W~ decay. Corrections to the 
LO results include O ^ ^ can be viewed as indicating 

the small size of part of the O corrections. Light fermion Yukawa couplings will be 

neglected throughout. 

The dehning assumption of universal theories translates into constraints imposed on the 
EFT. With proper held redehnitions, it is possible to transfer all the new physics effects into 
the bosonic sector. In other words, the effective operators generated by new physics involve 
SM boson helds only. Among all the bosonic operators respecting SM symmetries one can 
write down, the following ones contribute to —>■ W^W~ at LO: 


Ot = 

( 26 ) 

o„B = 

( 27 ) 

o^w = 

( 28 ) 

0,B = 

( 29 ) 

Ow = 

( 30 ) 

Ob = 

( 31 ) 

Ohw = tg{D^Hya^{D^H)W^,, 

( 32 ) 


( 33 ) 

6 ^ ^ 

( 34 ) 


We have adopted the notations in Refs. uniiig, and follow the conventions of Peskin and 
Schroeder [66] where — igW^^ — ig'B^Y, IE“^ = 9^IE“ — 

hP^ = CgZ^ -|- SgAfj^, Bfj, = —SqZ^ -|- CqA^j,. 

While our calculation of involves only the operators listed above, we emphasize that 

their Wilson coefficients do not have unambiguous meanings unless the full set of dimension- 
six operators chosen for the analysis is specihed. Usually a complete operator basis is 
desirable, although in other cases it is helpful to work with an over-complete operator set. 
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as long as one is careful about the RG running of the Wilson coefficients. A complete basis 
does not have to contain all the operators listed above. For example, the “EGGM basis” in 
Ref. [12] eliminates Ohw and Ohb via integration-by-parts in favor of two operators that 
affect Higgs physics only: 


Ohw —t Ow — ^(Oww + Owb) 

where 

Oww = 

(35) 

Ohb —t Ob — -^{Obb + Owb) 

where 

Obb=9'‘'\H\^B^,B^\ 

(36) 


Other operator bases, such as the one in Ref. na, trade some of the bosonic operators for 
those involving fermions via held redehnitions (or equivalently, equations of motion), which 
makes the dehnition of universal theories less transparent. We will keep all the operators 
above in the calculation, so that the hnal results can be easily adapted to and interpreted 
in different bases. 


B. Calculation of new physics effects 


There are hve aspects of the SM calculation of e’''e —>■ W^W polarized cross sections 
that are affected at LO by the operators listed above. 


• Modihcations of VWW vertices {V = 2’, 7 ). Assuming on-shell W^W , the vertex 
functions, dehned by the Feynman rules 


Za(9) ^ (37) 

Mq) ^ W;{k+)W-{k_) : zeFf "(g, k+, k_), (38) 


can be parametrized by g^-dependent form factors [39] 

(39) 

There are additional form factors if C and/or P violation is allowed, but these do not 
interfere with the LO SM contribution and are thus not considered here. In the EFT, 
the form factors read 


/;■(«=) - 1 + A7(,=) - 1 + Asr + ^ 


w 


+ 1 C2W, 


(40) 
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f^{q^) = Af^iq^) = Xv, (41) 

fs = 2 + A/3(g^) = 2 + Ag^ + Any H H —2 (^ (42) 

V Sg \zmyy / 


where 


AgJ = 

Agf = 

Ak^ = 

Akz = 

— 


0 , 


Ac^ 9 ^ 


{chW + Cw), 


--r^{cHw + Chb — ^Cwb) 
4s0 


A^f - -^Ak,, 


Az — 


-Czw 


( 43 ) 

( 44 ) 

( 45 ) 

( 46 ) 

( 47 ) 


are the commonly-used anomalous TGC parameters. We have rescaled the gauge 
helds hh“, to have their kinetic terms canonically normalized, and simultaneously 
rescaled g, g' so that gW^^, g'B^^ (and hence the gauge interactions of the fermions) are 
unchangec^ The weak mixing angle is redehned accordingly to retain sq = , ^ , 

y V 9 ^+ 9 '^ 

and the mass eigenstate fields are still dehned by = ceW^—sgBg^, = sgW^+ceBg^. 
The term , on the other hand, is not rotated away, which leads to noncanonical 

normalizations for Z„, as well as kinetic mixing. 


• Corrections to the s-channel Z /7 propagators. These can be viewed as corrections to 
the external leg of the VWW vertices: 

Tf ^ ^ ( 1 + , , Azziq^)) T^z^ + -A^ziq‘')Tf;’^\ (48) 

V q 

r7 ^ (1 + r7 + ^ (49) 

bQ q fiL^ 

where 

Azz{,q^) = ^ =- §-ct + + cw + cb) — %^{clc2w + sIc2b), (50) 

q‘^ q‘^ 2 

A'yziq^) = ^ 2^ ^ = Z - -^^{‘^CWB + CVK + Cb) — ^CgSg{c2W ~ C2 b), (51) 

q^ AcgSg 


This is possible at the dimension-six level in the EFT because the kinetic terms for and Wj^ are rescaled 
by the same factor. In other words, dimension-six operators do not generate a nonzero U parameter. 
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A^^[q ) _ 




= —^{4:CwB + Cw + Cb) - o{^0C2W + C0C2b) 


9 o \-^^yy^ ' ^yy o' 

1 


(52) 


are self-energy corrections due to new physics. Here and in the following, livv repre¬ 
sent the new physics contribution to the self-energies. Since in the SM, 


Eqs. (48) and (49) are equivalent to an additional contribution to the form factors: 

^2 


cZf2\ ^ A ^„2n , 


Ase/f (g^) = 


— m| 


^zziq ) H- \z{q ), 


ce 


nfr.A - A r^2\ , Co q ^ ,^2\ 




So q^ - m| 


^'yziq ), 


A../t (9^) = 0 . 

A.,/3^(,") = 2A../p4r"), 


(53) 

(54) 
(56) 
(56) 


where the subscript “se” stands for “self-energy”. 

Shifts in mw which enters the kinematics. With the Higgs VEV rescaled to leave 
the W boson mass term hE^hE“^ unchanged, only O 2 W contributes to mw via 
self-energy corrections 

^ (57) 


mfy —)■ -I- Uwwi'cnw) = '<ciw ( 1 ~ , 


0 

where mw should be understood as a shorthand for We emphasize that this step 
is essential regardless of whether mw is in the input observables set because the direct 
new physics corrections ^ are dehned with respect to the Lagrangian parameters 
(there are cancellations between direct and indirect contributions if mw is an input 
observable). There are also shifts in m^, but these are already contained in the 
propagator corrections calculated above. 

The held strength renormalization factors. The cross sections are simply rescaled. 


o'L,R —^ R 


(1 + — (Tl,r ~ ^*"2146^ • 


(58) 


Indirect contributions via shifts in the input observables. Using the input observables 
set we have 


1 

imz = + ^( 4 cwb + Cw + Cb) - ^^ 2^2 


{clc2w + sIc2b), 


(59) 


e-^e 
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= —^{^CwB + Cw + Cb)- 

^Gp = 0 because v has been rescaled such that YlwwiQ) = 0. 
Assembling all the pieces, we arrive at the dual result 


(60) 


0 crL,R — 


<^L,R 


da 


da 




iy 


df. 


dm^w 


(61) 


+ 2 n(^^y( 

Explicit expressions for aL,R and can be found in Appendix A Up to O cor¬ 

rections, we can trade the parameters in the hnal result of the calculation for the input 
observables {mz, Gp, cn} through LO relations, and use the LO results for (see Fig. 

for the size of NLO corrections to d^pp^i')- Numerically, using the reference values listed in 


Section III for the input observables, we obtain, for s/s = 200 GeV, 


l = 0.0192[chw + cw ) + 0 . 00345 ( ci^vv + chb ~ '^cwb) + 0 . 00667 c 3 vi / 

—0.0967(4cvi/i3 + Cw + cb) + 1.66ct — 0.1S3c2w + 0.0442c2b, (62) 

= —1.32(cbw + ^w) + 0.640(c/fvv + cpB ~ ^cwb) ~ 0.0898c3w 
-|-1.67(4cwi3 + Cw + Cb) — 5.45ct — 0.17302^ — l-49c2s, (63) 

= 0.00513(c/f^y + Cw) + 0 . 0101 (ci^vv + chb ~ 4cvKij) + 0.00566c3vi/ 


—0.0782(4cw_b a Cw + cb) + 1.58cj' — 0.183c2w A 0.0281c2i3, (64) 

and for y/s = 500 GeV, 

= 0 . 0835 ( cj^vp + cw) + 0.0277(chw + chb — dciys ) + 0 . 0191 c 3 ^y 

—0.115(4ciyi3 + Cw A Cb) + 2.3Sct — 0.253c2w + 0.0396 c2_b, (65) 

d^^aji = —3.2b{cHw + cw) + Q.Q^{chw + chb ~ 4cwi?) ~ 0.0426c3^ 

+8.41(4cvi/_b + Cw + Cb) — 2.61ct' — 0.0826c2w — 8.33 c2_b, (66) 

= 0.0497{cfjw 4“ Cw) + 0.0546(^chw 4- chb — 4lCwb) 4- 0.0189c3iy 
—0.0804(4c^^^ A Cw A cb) A 2.3Gc'b — 0.252c2w A 0.00563c2ij. (67) 


We have also shown the new physics contributions to the unpolarized cross sections in the 
equations above, which are directly calculated from 

( 68 ) 

o'l 4- aji o-L A aji 
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The results for the left-right asymmetries are derived from 


XNP /t _ 

0 ^LR — 0 (^R — 0 cr, 


NP „ xNP, 




-0.0214 for = 200 GeV 


1 - A 


LR 


(69) 
(LO), (70) 


LR 


-0.00818 for = 500 GeV 


and will not be listed explicitly. 


C. Interpretation of resnlts 


To interpret these results, namely to see the role played by e’''e“ —)■ W^W~ in probing 
the EFT parameter space, we need to go to a specihc basis, and compare our results with 
other experimental constraints on the Wilson coefficients in this basis. We consider the 
EGGM basis [I2] for illustration, and relate four combinations of Wilson coefficients in this 
basis to the oblique parameters S, T, W, Y [6T1 [671 168] : 

(71) 

(72) 

(73) 

(74) 


A aS 00.^, . 

S = — —— 113 ^( 0 ) — ^(4cr/b + cw + cb), 


4s^ 
f = aT = 


S0 


1 


m 


w 


[niyv(/(0) — 1133(0)] — Cr, 


^ 33 ( 0 ) - 772'^2W, 
e 




2 - br (0)-4^2^25- 


Y = 


We have adopted the definitions in Ref. [TB], which differ from Ref. [HI] by sign. In universal 
theories, these four parameters are sufficient to describe the Z-pole data at LEPl, mw 
measurements, and e'''e“ —)■ // data at LEP 2 | 6 T]. Each of them is constrained at the 
10“^ level [61]. It is clear that these oblique parameters are not sensitive to C 3 vl- Also, they 
constrain only one combination of the three coefficients cwB,cwyCB- Ghoosing the other 


two combinations to be cwb and cw — cb, we rewrite Eqs. (62 67) in the EGGM basis as 
follows. 


^^^crr(v/s = 200 GeV) = -0.0521cwb + 0.00958(cw - cb) + 0.00667c3w 

-0.826A + 1.66f-1.74IV +0.419V, (75) 

^^^cpr(v/s = 200 GeV) = 0.0805cvkb - 0.660(ciy - cb) - 0.0898c3w 
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+9.53^ - 5.45f - IMW - U.IY, (76) 

= 200 GeV) = -0.0507cvKij + 0.00256 (cm^ - cb) + 0.00566c3vf 

-0.7175 + l.SSf - 1.74PG + 0.267F, (77) 

(5^^t^L(ys = 500 GeV) = -0.278 ch.s + 0.0418(cvi/ - cb) + 0.0191c3vk 

-0.6955 + 2.38f - 2.401V + 0.375V, (78) 

5'^^cTi^(A/s = 500 GeV) = —IO.Icvkb ~ 4.12(cpi/ — cb) — 0.0426c3vi/ 

+40.75 - 2.6lf - 0.7831V - 79.0V, (79) 

(5^Pa(ys = 500 GeV) = -0.318cvt^B + 0.0249(ch/ - cb) + 0.0189c3vi/ 

-0.5275 + 2.36f - 2.391V + 0.0534V. (80) 


It should be noted that the coefficients of 5 in these equations are not unique, as they depend 
on the choice for the other two combinations among cwb,cw,cb- Our choice is motivated by 
the observation that in weakly-coupled new physics scenarios, and Ob are “potential- 
tree-generated”, while OivB is “loop-generated” [TUI[69H73] . In the limit \cwb\ "C |cvy|, |cb|, 
the 5 parameter involves cw + cs, while the orthogonal combination cw — cb can be probed 
by e+e- ^ IV+IV. 


Eqs. (75 80) demonstrate the complementarity of e'''e 1V’''1V and other precision 


data in probing the parameter space of universal theories at the observables level. As such 
they are free from ambiguities and additional assumptions. It is seen that cross sections 
at \/s = 500 GeV tend to have stronger dependence on the Wilson coefficients, and hence 
better sensitivity to new physics effects. This is not surprising since there are contributions 
that scale as a fact that is easier to see in a basis where operators involving fermions 
are retained |6^. We have chosen to focus on bases that maximize the use of bosonic 
operators for easier connection with the literature on universal theories, but one should not 


misinterpret the relatively strong dependence on multiple Wilson coefficients in Eqs. (75- 


80) as necessarily indicating experimental sensitivity to all of them. To draw quantitative 


conclusions about the constraining power of —)■ 1V+1V“, a global analysis should be 

carried out, where correlations among the constraints on different Wilson coefficients become 
clear. The calculations and results presented in this paper can be used as a starting point. 
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V. CONCLUSIONS 

Precision electroweak analyses will continue to contribute to our understanding of Nature 
in the post-Higgs-boson-discovery era. A truly consistent approach to precision analyses, 
both as a consistency test of the SM and as an indirect probe of new physics, requires not 
only the experimental and theoretical precision, but also the dependence of the calculation on 
the input observables to be well understood. Though in many cases parametric dependence 
is insignihcant compared with experimental and theoretical uncertainties, it should not be 
taken for granted that it will remain so, especially when the projected experimental precisions 
greatly exceed the current ones. Instead, careful justihcation is needed. 

In this article, we have investigated the role of —)■ W~^W~ in the precision program, 

motivated by the projected per-mil-level cross section measurements and the complemen¬ 
tarity between this and other processes in probing new physics effects. The latter is also 
relevant for precision studies of the Higgs boson. We have utilized the best SM calculations 
available as public codes, and presented the results for several inclusive observables in terms 
of expansion formulas. From these one can directly read off the dependence on the input 
observables, which allows us to justify the exclusion of inclusive e+e“ —)■ W~^W~ observ¬ 
ables from the global precision electroweak analysis as a test of the SM. In fact, even in the 
future, the experimental uncertainties are still expected to dominate over the parametric 
uncertainties, indicating an insignihcant contribution to the function. Our analysis also 
justihes the neglect of SM parametric uncertainties when using precision data to constrain 
new physics contributions above the 10“^ level. 

On the other hand, to demonstrate the interplay between e’''e“ —)■ W~^W~ observables 
and other precision data in probing new physics effects, we have considered a representative 
class of models, the universal theories. These models have been extensively studied in 
the oblique parameters framework, where four parameters suffice to incorporate the most 
stringent experimental constraints. For e+e“ —)■ W~^W~, we have gone beyond the TGC 
parametrization, and calculated well-dehned physical observables in the EFT framework, 
again using the expansion formalism. While some terms in the expansion can be mapped to 
the oblique parameters, other terms show the additional directions in the EFT parameter 
space accessible to —)■ W~^W~. Our calculations and results provide necessary tools 
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for a consistent global analysis, which will tell us where we are in the EFT parameter space 
if deviations from the SM are found in upcoming precision measurements, and will thus give 
us guidance for where to look for new physics if it fails to show up directly. 
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Appendix A: Leading-order results for e+e — )■ W~^W cross sections in the SM 


At LO, the full e+e“ —)■ 4/ process factorizes into production and decay. The 

decay branching fractions of W~^ (and same for W~) simply follow from hnal state count¬ 
ing, and are 1/3 for hadronic hnal states and 1/9 for leptonic hnal states. In the PDG 
notation [BUI, the production cross sections for polarized e~ read 


f*t — 


( 2 2 2 2 \ 
aL(s;m^,m^y,e ,Sq) = 


Stts^ 


dt 


't+ 


2sl-l 


2sl s — 




2sl-l 


2sg s — 


0 -^( 5 ; m|,m^,e^) = 


Stts^ 


df 


'*+ 


rur. 


1 ) + ■^E{s,t;rnw) 


-1 A(s,t;m^), 


,(A1) 

(A 2 ) 


where t± = 


ill 


1 ± 


Ami, 


— 1. The functions A, /, E come from the s-channel, 
interference, and f-channel contributions, respectively. Explicitly, 


Af 2 \ s'^ f t S 




E{s,t;m‘^w) = 



Vi 

1 I' 




4m^ 


, ( 

'As 

+ 2s 1 


S2 

ft. 


6t 
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3m^ 


2 m^ 
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w 


4m^ s 


st ’ 
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2rn^ 


3t 2t^\ 1 / 20t 12^2 

2+- + ^ --17+ — + — 
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2 H- \- — 
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5 
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2 + - 
s 




(A3) 


2m^ 
t E 


(A4) 


•(AS) 
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Eqs. (A1) and (A2) are written in terms of the parameters that directly enter the calcnla- 
tion, all of which are not independent. This parametrization is convenient for the calcnlation 
of direct contribntions from new physics, namely ^aLR- Bnt regarding the indirect contri- 
bntions — 'Yhv R,i'ii', one shonld be carefnl not to over-connt the nnmber of independent 
inpnts. Using {mz,GF,o:} as inpnt observables, we have 


where 


2m| / dai dai drri^ dai dsj \ 
ai v9m| dm^ dm\ ds^ dm^J ’ 
Gf f d(TL dml^r d(TL dsl \ 
ctl \dm\r OGf dsl ^Gf) ’ 

2 a / dai dmly dai dsl \ 
ai \dm^ da dsl / ’ 


Gfi a) 


sl{m\,GF,a) 


/i ^ 

2 + Y ^2GFml J ’ 


1 

2 



' 1 - 


47ra 


V^GFm^^ 


(A 6 ) 

(A7) 

(AS) 

(A9) 

(AlO) 


are nsed to calcnlate the partial derivatives involved. For ctr, jnst drop the terms with si- 
It is nnderstood that any dependence on the calculational inpnts in the hnal resnlts shonld 
be traded for the three inpnt observables, which can be done at LO, before nnmerical valnes 
for the inpnt observables are plngged in. These LO resnlts for the expansion coefficients are 
also used for comparison with the results from RacoonWW in Fig. 


It is shown in Sec. IV that a large part of new physics effects is effectively characterized 


by the form factors dehned in Fq. (39). The partial derivatives in Fq. (61) read 
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where ^ are obtained by setting = fi and taking derivatives, 
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Appendix B: Technical details of parametric dependence calculations 


Among the hve observables introdnced in Sec. [TT| only two are independent, which we 
choose to be a and an for the present calcnlation. The reference valnes and expansion 
coefficients for the other three observables can be derived from those for a, as follows: 


^ref 

= 

O^ref ^ref 

2cr - cr^ , 

CaL,i' 

= Ar(2^’ 
af 

4 ref _ 
^LR — 

ref 


. 

^ref’ 

Ulr 





-rpf 

4 — 

^LR “ 

^ref’ 


= Co-R,i' ~ 


ref _ ref \ 

^a,i' ^R )'! 

(Co-,j' — Co-^^j/), 


(Bl) 

(B2) 

(B3) 


We nse the RacoonWW package to calcnlate e+e“ —>■ ud^~9^ cross sections a, a^, 
adopting the recommended choices for the switches in the inpnt file bnt removing all sep¬ 
aration cnts. The program is rnn mnltiple times with shifted inpnt observables before the 
expansion coefficients Ci^i' are calculated by the hnite difference method. 
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+ O [f] 
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O, 

^i,mz 


hi^myy 

h- 

hi,Os 

hi^rnjj 

a 

0.03 

0.03 

0.03 

0.1 

0.1 

0.1 

<yR 

0.005 

0.03 

0.005 

0.1 

0.1 

0.05/0.1 


TABLE IV. The choices for hj j/ in Eq. (B4), the fractional shift in the input 0-, with respect to 
for the extraction of For 0-0^ is used for = 200 GeV while 0.1 is used for 

a/s = 500 GeV. 


This is clearly CPU-time-consuming since the calculation involves Monte Carlo integra¬ 
tion over the final-state four-fermion phase space. We thus keep the number of weighted 
events generated per run of the program iVevents at a minimum that still allows interesting 
information to be extracted. In particular, for the extraction of com¬ 

bine the results from the two branches “slicing” and “subtraction” of Racoon WW, with 
^events = 10^ for each branch, which is the recommended minimum. The two branches differ 
in the treatment of IR singularities, and they are found to yield compatible results. For 
the much smaller ^^t Nevents = 5 x 10^ to further reduce Monte Carlo 

error, while running the faster “slicing” branch only. The reference values are also 

obtained with the “slicing” branch with iV/vents = 5 x 10^. 


The choices of in Eq. (B4) are listed in Table IV These h’s are small enough so 
that reducing h further leads to larger Monte Carlo error bars that encompass the current 
ones. In other words, the truncation errors due to the discarded O (h^) term are negligible 
compared with the Monte Carlo errors. Of course, the optimal h’s would be smaller if the 
Monte Carlo errors are reduced by increasing A/vents- The Monte Carlo errors quoted in 
Tables and are obtained assuming individual runs of RacoonWW yield uncorrelated 
errors. This assumption is probably very conservative because, since the same random 
number seed is used for each run, the results for a, ur for slightly different values of the 
input observables tend to fluctuate in a correlated way. Therefore, the central values of the 
expansion coefficients shown in Tables [I] and [IT] are likely to be more robust than the error 
bars suggest. 
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